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Abstract 

The structure of the solutions to the one dimensional focusing nonlin¬ 
ear Schrodinger equation (NLS) for the order N in terms of quasi rational 
functions is given here. We first give the proof that the solutions can be 
expressed as a ratio of two wronskians of order 2N and then two determi¬ 
nants by an exponential depending on t with 2N — 2 parameters. It also 
is proved that for the order N, the solutions can be written as the prod¬ 
uct of an exponential depending on t by a quotient of two polynomials of 
degree N{N -I- 1) in a; and t. The solutions depend on 2N — 2 parame¬ 
ters and give when all these parameters are equal to 0, the analogue of 
the famous Peregrine breather Pn- It is fundamental to note that in this 
representation at order N, all these solutions can be seen as deformations 
with 2N — 2 parameters of the famous Peregrine breather Pn- With this 
method, we already built Peregrine breathers until order N = 10, and 
their deformations depending on 2N — 2 parameters. 


1 Introduction 

The term of rogue wave was introduced in the scientific community by Draper 
in 1964 [1]. The usual criteria for rogue waves in the ocean, is that the vertical 
distance from trough to crest is two or more times greater than the average 
wave height among one third of the highest waves in a time series (10 to 30 
min). The first rogue wave recorded by scientific measurement in North Sea 
was made on the oil platform of Draupner in 1995, located between Norway and 
Scotland. Rogue waves in the ocean have led to many marine catastrophes; it 
is one of the reasons why these rogue waves turn out to be so important for the 
scientific community. It becomes a challenge to get a better understanding of 
their mechanisms of formation. 

The rogue waves phenomenon currently exceed the strict framework of the study 
of ocean’s waves and play a significant role in other fields; in nonlinear optics 
[2], Bose-Einstein condensate [3], atmosphere [4] and even finance [5]. 

Here, we consider the one dimensional focusing nonlinear Schrodinger equation 
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(NLS) to describe the phenomena of rogue waves. The first results concerning 
the NLS equation date back the works of Zakharov and Shabat in 1968 who 
solved it using the inverse scattering method HIT]. The case of periodic and 
almost periodic algebro-geometric solutions to the focusing NLS equation were 
first constructed in 1976 by Its and Kotlyarov 011]. In 1977 Kuznetsov found 
the first breather type solution of the NLS equation m-, a simular result was 
given by Ma m in 1979. The first quasi rational solutions to NLS equation 
were constructed in 1983 by Peregrine [T^. In 1986 Akhmediev, Eleonski and 
Kulagin obtained the two-phase almost periodic solution to the NLS equation 
and obtained the first higher order analogue of the Peregrine breather [T^ . 
Other analogues of Peregrine breathers of order 3 were constructed and initial 
data corresponding to orders 4 and 5 were described in a series of articles by 
Akhmediev et ah, in particular in |141115] using Darboux transformations. 
Quite recently, many works about NLS equation have been published using dif¬ 
ferent methods. In 2010, rational solutions to the NLS equation were written 
as a quotient of two wronskians m- In 2011, the present author constructed 
in m another representation of the solutions to the NLS equation in terms of 
a ratio of two wronskians of even order 2N composed of elementary functions 
using truncated Riemann theta functions depending on two parameters; ratio¬ 
nal solutions were obtained when some parameter tended to 0. In 2012, Guo, 
Ling and Liu found another representation of the solutions as a ratio of two 
determinants m using generalized Darboux transform; a new approach was 
proposed by Ohta and Yang in [20] using Hirota bilinear method; finally, the 
present author has obtained rational solutions in terms of determinants which 
do not involve limits in |21] depending on two parameters. 

With this extended method, we present multi-parametric families of quasi ratio¬ 
nal solutions to the focusing NLS equation of order N in terms of determinants 
(determinants of order 2N) dependent on 2N — 2 real parameters. With this 
representation, at the same time the well-known ring structure, but also the 
triangular shapes also found by Ohta and Yang j^O], Akhmediev et al. [IS] are 
given. 

The aim of this paper is to prove the representation of the solutions to the 
focusing NLS equation depending this time on 2N — 2 parameters; the proof 
presented in this paper with 2N — 2 parameters has been never published. This 
is the first task of the paper; then we deduce its particular degenerate represen¬ 
tations in terms of a ratio of two determinants of order 2N. The second task 
of the paper is to give the proof of the structure of the solution at the order N 
as the ratio of two polynomials of order N{N -1-1) in cc and t by an exponen¬ 
tial depending on t. This representation makes possible to get all the possible 
patterns for the solutions to the NLS equation. It is important to stress that 
contrary to other methods, these solutions depending on 2N — 2 parameters give 
the Peregrine breather as particular case when all the parameters are equal to 0 
: for this reason, these solutions will be called 2N — 2 parameters deformations 
of the Peregrine of order N. 

The paper is organized as follows. First of all, we express the solutions of the 
NLS equation using Fredholm determinants from these expressed in terms of 
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truncated functions theta of Riemann first obtained by Its, Rybin and Salle 
[9]; the representation given in theorem 2.1 is different from those given in [^. 
From that, we prove the representation of the solutions of the NLS equation in 
terms of wronskians depending on 2N — 2 parameters. We deduce a degenerate 
representation of solutions to the NLS equation depending a priori on 2N — 2 
parameters at the order N. 

Then we prove a theorem which states the structure of the quasi-rational solu¬ 
tions to the NLS equation. It was only conjectured in preceding works [13 [H]. 
Families depending on 2N — 2 parameters for the N-th order as a ratio of two 
polynomials of x and t multiplied by an exponential depending on t are obtained; 
it is proved that each of these polynomials have a degree equal to N{N + 1). 


2 Expression of solutions to the NLS equation 
in terms of wronskians 


2.1 Solutions to the NLS equation in terms of 9 functions 


For r = 1,3, we define 


dAx,t) = Efcetopp'v exp 


+ 2d^t + 


In 


7 .^- 7 ^ 


-I- Trie„ + e, 


0 J” ’ 


( 1 ) 


In this formula, the symbol denotes summation over all 2N-dimensional 

vectors k whose coordinates ki, are either 0 or 1. 

The terms Kjy,(5i/,7jy and Xr^u are functions of the parameters Ajy, 1 < v < 2N; 
they are dehned by the formulas : 


= (r - l)ln^|^, r = l,3. 


The parameters —1 < Aj/ < 1, = 1,..., 27V, are real numbers such that 


( 2 ) 


— 1 < Xn+1 < Xn+2 < ■ ■ ■ < X 2 N < 0 < Atv < Ajv-i < ... < Ai < 1 , . 

Xn+j=-X„ j = 

The condition (|31) implies that 

^j+N — ^j+N — 7j 1 — Xr^j-, j — 1; • ■ ■ ; Xk. (4) 

Complex numbers e,y 1 < < 2N are dehned in the following way : 

ej = ittj — bj, CN+j = iaj +bj, 1 < j < TV, a, 6 S R. (5) 

e,y € {0; 1}, (p, z/ = 1... 2TV are arbitrary real numbers. 

With these notations, the solution of the NLS equation 

ivt+vxx + ‘2.\v\‘^v = 0, ( 6 ) 
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can be expressed as m) 


vix, t) = exp(2ii - i^p), 


(7) 


2.2 Prom 6 functions to Fredholm determinants 


To get Fredholm determinants, we have to express the functions 9r defined in 
CD in terms of subsets of [1, 2N] 


9r{x, t) 


E IK-D- If 


Iv - Jfj. 


X 



— 2S,yt + Xr,v 



( 8 ) 


In ([8]), the symbol X]jc{i 2 N} denotes summation over all subsets J of indices 
of the set {1,2iV}. 

Let I be the unit matrix and Cr = {cjk)i<j,k< 2 N the matrix defined by : 

/ Y\r]pfi Tyryl + i I \ 

= (-1) TT I -r exp(zK,.a; - 25^1 + Xr,v, + e^), (9) 

\Ar,Pu 


ej=j l<j<N, e,=N + j, N+l<j<2N. (10) 

Then det(/ + Cr) has the following form 


det(/ + Cr) 


n iK-i)- n 

Jc{i,...,2W} J ugJ iJ.^J 


+ 7m 

Iv-lki 


exp(zKi^x 


25i,t + Xr^u + Cl/). (11) 


Comparing this last expression CID with the formula ([8]) at the beginning of 
this section, we have clearly the identity 


9r = det(/ + Cr)- 


( 12 ) 


We can give another representation of the solutions to NLS equation. To do 
this, let’s consider the matrix Dr = {djk)i<j,k< 2 N defined by : 


dufi 


(- 1 )" n 


7r; -7 m 


exp(zKi,a: 


2Sjrt Xr^iy “t“ (-r/). 


(13) 


We have the equality det(/ + Dr) = det(J + Cr), and so the solution of NLS 
equation takes the form 


v{x,t) 


det(/ + Dz{x,t)) 
det(/ + Di{x, t)) 


exp(2zt — ip). 


(14) 
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Theorem 2.1 The function v defined by 


v{x, t) 


det(/ + Dz{xfi)) 
det(/ + Di{x, t)) 


exp(2i< — iip). 


(15) 


is a solution of the focusing NTS equation with the matrix Dr = {djk)i<j,k< 2 N 
defined by 


du,={-If' n 


Iv-Tu 


exp{iKL,x 


‘2iSiit Xr^jy ef. 


where k^, 5^, Xr,v, Iv, Gu being defined m|0), Q) and (0). 


2.3 Prom Fredholm determinants to wronskians 

We want to express solutions to NLS equation in terms of wronskian determi¬ 
nants. For this, we need the following notations : 

= sin 0r,i/, i < V < N, 4>r,v = COS 0r,j/, iV -|- 1 < ^ < 2N, r = 1,3, (16) 

with the arguments 

0r,u = KvXf -I- iSiyt — iXr,vf + '^vV — i^vf, ^ < 2N. (17) 

We denote Wr{y) the wronskian of the functions f,!, ■ ■ ■, 4>r,2N defined by 


Wr{y) = det[{d!^ Vr,i/)i/,,ie[i,...,2Ar]]- (18) 

We consider the matrix Dr = (di//i)j/, ^e[i,..., 2 Ar] defined in m- Then we have 
the following statement 

Theorem 2.2 


det(/ -f Dr) = kr{0) X Wr{4>r,l, ■ • ■ , fir,2N){0), 


(19) 


where 


kr{y) 


2^^ exp{tj:l^, Qrf 

nl^2u;i\{i^-iuy 


Proof : We start to remove the factor {2i)~^e^^^’" in each row i/ in the 

wronskian Wr{y) for 1 <v < 2N. 

Then 


2N 

Wr = W e^^'-’f2i)-^{2)-^ X Wr, (20) 
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with 


(l_e-2*e.,i) + ... (i-^i)2W-i(l + (_l)2iVg-2*e.,i) 

(1 - 172(1 + ... ( z 72)^^"^(1 + 

Vhr = . 

(1 _ g-2ie,,2N) j^2Ar(l + e-2i©r.2iv) (•_;^)2ng-2ie,.2iv) 

The determinant Wr can be written as 


Wf — (X^ti^CXjfzCj ~h Pjk'jj 

where ajk = {-l)’"{ijj)’"~^, ej = and j3jk = (*7i)^”\ ^ < j < N, 

l<k< 2N, 

ajk = (-l)''“^(z 7 j)''“\ Cj = and Pjk = (*7i)''”\ + 1 < J < 27V, 

1 < fc < 27V. 

We want to calculate Wr- To do this, we use the following Lemma 

Lemma 2.1 Let A = (a^-)*_jg[i^,„_Ar], B = (%)*, je[i,...,Ar], 

(iVy )i jg[i jv]: the matrix formed by replacing in A the jth row of A by the ith 
row of B Then 


det{aijXi + bij) = det(aij) x (let{6ijXi 


det(iLij), 


Proof : We use the classical notations : A = (dji)i, je[i,....Az] the transposed 
matrix in cofactors of A. We have the well known formula A x A = det Ax I. 

So it is clear that det(24) = (det(24))'^“^. 

The general term of the product = {aijXi + ^ (hji)i 76 [i...,JV] 

can be written as 

Czj — “t“ bis^ X Ujs 

— 1 ^is^js T 1 bisOjs 

= Sij det(24)xi + det(iLy). 

We get 

det(cij) = det(aya;i + bij) x (det(24))^“^ = (det(24))^ x det(%a;i + 

Thus det(aijXi + bij) = det(24) x det{SijXi + )■ 

□ . 

We denote U = (aii)i,je[i,..., 2 Af]) V = (Ai)i,iG[i,..., 2 Af]- 
By applying the previous lemma, one obtains : 


Wr — det(aijei + /7y) 

= det(aij) X detiSyCi + ^1^) = det(U) x det(Syei + 

where (-ffyOi, je[i,...,iv] is the matrix formed by replacing in U the jth row of U 
by the ith row of V defined previously. 

The determinant of U of Vandermonde type is clearly equal to 

det(17) = i'ri-jm)- (23) 

2N>l>m>l 
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To calculate determinant Wr, we must compute now det(JIij). To do that, two 
cases must be studied : 

1. For 1 < j < JV. The matrix Hij is clearly of the VanderMonde type where 
the j-th row of ?7 in {7 is replaced by the i-th row of V. Clearly, we have : 

det(i7y) = (_l)Ar( 2 Af+l)+Af-l(i)Af( 2 iV-l) X M, (24) 

where M = M(mi,. .., rn 2 N) is the Vandermonde determinant defined by mk = 
7 fe for k ^ j and rrij = —yj. Thus we have : 

det{Hij) = —X 02Ar>i>fc>i, {mi - mk) 


— X Yl2N>l>m>lJ^j,m=/^j('yi- 7m) 7i 


t) X n/>i(7i+7*), (25) 


= (-l)J(i)^(2^-l) 


X Yl2N>l>m>l,l^j,m=iij('l'l 7m) X Y\i^j{"fl +7i)- 


To evaluate Wr, we must simplify the quotient qij := : 

_ (-W ni^j(7i+7i) _ (-W ni^j(7i+7d _ ni^j(7i+7i) 

“ ni<j( 7 j- 7 i)rii> 3 ~ (- 1 )^"^ 1115^3(71-73) “ 01^3(71-73) ■ 


(26) 


We can replace qij by defined by ~ > because det{SijXi + d^et^A)^ ) ~ 

det{dijXi + (similar matrices). 

We express rij in terms of absolute value; as j S [1;77] and 0 < 71 < ... < 
7 Ar < 1 < 72 Ar < ... < 7 Ar+i, we have : 


= 11 / 3 ^ 3 ( 7 /+ 7i) = It+ 7iI• (27) 


So the term Vij can be written as 


_ / 3 \3 1115113 l71“t"7i| 

'' ni5^ii7i-7ii 


r 


_.| Yd) ni5^3lT'i+7i| ^ ^ .g-2i0,,,i(O) 


(- 1 ) 


Ill5^il7i-7i| 


with respect to the notations given in (fTUl) and (|T^ . 


(28) 


2. The same estimations ior N + 1 < j < 2N are made; det Hij is first 

det(i733) = (_i)^(2Af+i)+Af-i(i)Af(2Ar-i) x M, (29) 

with M = M{mi ,..., m 2 N) the Vandermonde determinant defined by mu = 7 fc 
for k ^ j and mj = — 73 . Thus we have : 

det(i7y) = (i)^( 2 ^-i) X n 2 Ar>/>fe>i. iw-i -™fc) 

= (z)^(2iV-l) X n 2 Ar>J>m> 17 / 3 .m/ 3 ( 7 / -7m) X n/<i(-7i “T) X n;>3(T +T), (30) 
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The quotient := equals : 


qij = 


(-i)"-i(i)"(2"-i>xn 


2N>l>m>l, l^j, Tn5^j(^'< lrn)'xY\l^jill'^li) 


ni<j(7j-7i) ni>i(7j-7j) (-i)j-i 1115 ^^( 71 - 73 ) ni^j(7i-73)' 


(31) 


for the same reason as previously 


We replace qij by r„- defined by 
exposed. 

rij is expressed in terms of absolute value; as j G [N + \]2N] and 0 < 71 < 

... < 7 iv < 1 < 72 Ar < • ■ • < 77 V+ 1 , we have : 

-7i) = (-1)2^-*+^It -lA, Yh^jin + 7 *) = Oz/j It + 7*1 • (32) 

So the term rij can be written as 


^ _ { _ 1 \N-\-i ritT^j I _ /_ 1 '\e(z} l7t+7^| _ — 22 ©^ i(0) 

T ni^il7i-7i| “ n,^il7i-7d “ 

with respect to the notations given in (fTOll and (IT^ . 

Replacing a by det 114 can be expressed as 

detll4 = det(17) x det{Sijei + “ det(17) x det(^yei + rij) 


(33) 


= det(t/) n?=i det(% + (-l)dd 


det(i7) 

1-1 Id*) ■ 


7^+7^ 

ii-'yi 


e2*e.,.)_ 


(34) 


We estimate the two members of the last relation (IM)l in y = 0, and using 
we obtain the following result 

det W 70 ) = n 2 ZV>z>***>i(T - 7m) O-fi 


xdet(% + (-l)ddn,^, S 

= iN(2N-l) q2Ar^ n^-7^, _ 


„ 2 *e. 


i(0)) 


(35) 


' n ■:2 nS(T - T)e-^*^- det(% + c,,) 

= l^(2zv-i) UiZihj - 7*)e-"*^'=i®-do) det(/ + a) 

= z^(2zv-i) nCi'(7, - 7*)e-"*^" ®-do) det(/ + D^). 

Therefore, the wronskian 114 given by (EOl) can be written as 
Wricl}r,l,.-.,Ar,2N){0) = n'fl 6*®-(0) (2)-21V (^)x ITr 

= n-fi e*©-d0)(2)-2W(,)-zv^zv(27v-i) det(/ + D^) (36) 

= (2)-^^ n-^2 nCi'(T - T)e-^^- det(/ + D^). 

As a consequence 


det(/ + Dr) = kr{0)Wr{(j)l,. ■. ,(/'2Zv)(0). 


(37) 

















2.4 Wronskian representation of solutions to the NLS equa¬ 
tion 


From the initial formulation m we have 


v{x, t) 


det(/ + D^(x,t)) 
det(/ + Di{x, t)) 


exp(2jt — lip). 


Using (HU), the following relation between Fredholm determinants and wron- 
skians is obtained 


det(/ + D3) = ^3(0) X ■ ■ ■ , (/)r,2Af)(0) 


and 

det(/ + Ds) = A:3(0) X lU3((/)r,l, . . . , (/>r,2Ar)(0). 

As 03^j(O) contains N terms x^^j 1 < j < and N terms —Xsj 1 < j < A^, we 
have the equality ^3(0) = fci(O), and we get the following result : 


Theorem 2.3 The function v defined by 


v{x,t) 


1U3(^3,1, ■ ■ ■; </>3,2jv)(0) 
Wi{(j)ip, . . . , fil.2N){0) 


exp(2if — i(p). 


is a solution of the focusing NLS equation depending on two real parameters a 
and b with fil, defined in m 


4>r,u = siti(kvx/2 + i5vt — iXr^vf^. + juV — i&vl2L)^ 1 < < A^, 

(f>r,v = cos(Kjya;/2 + iS^t — ixr^vj^ + — ie^j^), N + 1 < v < 2N, r = 1,3, 

Kj/, (5^, Xr^i^, Cu being defined in^, (0) and 


3 Families of multi-parametric solutions to the 
NLS equation in terms of a ratio of two de¬ 
terminants 

Solutions to the NLS equation as a quotient of two determinants are constructed. 
Similar functions defined in a preceding work m are used, but modified as ex¬ 
plained in the following. The following notations are needed : 

Njy — KuX j2 “t“ i5i/t ix^^u f 2 %ei/l2^ 

Yl, = k^xI2 -f — ix\^vl^ — ie-vj^i 

for 1 < 13 < 2N^ with 8^, Xr,v dehned in ©■ 

Parameters Cy are defined by ®. 

Here, is the crucial point : we choose the parameters and bj in the form 

N-l N-l 

aj = Y, bj = Y 1<J<N. (38) 
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Below the following functions are used : 


(P4j+l,k = lt smXk, (f4j+2,k = it cosXk, 

4j + l ■ 4j+2 

^4j+lk = -Ik SmXk, (/?4j+4.fc = -7fc COSAfc, 

for 1 < fc < -/V, and 

2N-ij-2 V 2Af-4i-3 • ^ 

V54i+i,Af+fe = 7fc cosA^r+fc, i^4i+2,A+fe = ~7fc sinA^r+fc, 

2Ar-4i-4 2A-41-5 • ^ 

'^4i+3,A+fe = ~lk COSAjv+fc, i^4i+4,A+fe = Ik SinA^r+fc, 

for 1 < fc < TV. 

We define the functions ipj^k for 1 < j < 2A^, 1 < fc < 2N in the same way, the 
term Xk is only replaced by Yfc. 

■04i+i.fe = lt^~^ sin Yfc, ^4i+2,fc = it^ cosYfe, 
i’4j+3,k = sin Yfc, 'Ip4j+4,k = -it^^"^ cos Yfc, 

for I < k < N, and 

, 2N-4j-2 , 2N-4j-3 ■ 

■(/’4i+l.A+fe = 7fc COsYAT+fc, ■i/'4i+2.A+fe = “7fc SluYAr+fc, , . 

, 2A-4j-4 -5^ , 2A-4j-5 • 

V’4i+3.A+fe = “7fc COsYAT+fc, V’4i+4.A+fe = 7fc SluYAf+fc, 


ioT 1 < k < N. 

Then it is clear that 


can be written as 


q{x,t) := 


1Y3(0) 

fYi(O) 


<,(*4) = ^ = (43) 

Ai det('iyj,fc)j_fcg[i_2Af] 

We recall that Xj = 1 —2je^. All the functions (pj^k and i/'i.fc and their derivatives 
depend on e and can all be prolonged by continuity when e = 0. 

Then the following expansions are used 

V3yfc(x,t,e) = + (x, t,0), 


i[0] = V3yi(x,t,0), l<j<2iV, l<fc<fV, l</<iV-l, 

Pj,N+kix,t,e) = Yt, 7^‘/5i.A+iWfc^'e^'+0(e2^), Pj,N+i[l] = —|^^^(x,t,0), 

(/?y7v+i[0] = V3j,A+i(x,t,0), 1<J<27V, l<k<N, 1 < I < N - 1. 

We have the same expansions for the functions 'fpj^k- 

1 3‘^^'ih ■ 

i;j,k{x,t,e) = Y + V’j.iW = {x,t,0), 
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l<j<2iV, l<k<N, 1<1<N-1, 


1pj,N+k{x,t,€) 


^ ^3,N+i[1] 

1^0 ^ '' 


3^21 


{x,t,0), 


ipj^i^^i[0] = tj;jj^^i{x,t,0), i<j<‘2N, 1 < fc <-/V, A^ + l<fc<2A^.. 

Then we get the following result : 


Theorem 3.1 The function v defined by 


v{x, t) = exp(2it 


i(p) X 


det((njfe)^ ) 

det((djfc)3_fc(i[i^2N-]) 


(44) 


is a quasi-rational solution of the NLS equation (0) 


ivt + Vxx + 2|w|^z; = 0, 


where 

o2fc-2 

0), 1 < j < 2iV Ujk = Q^ 2 k- 2 ^ {x,t,Q), 

a2fc-2 

njN+i = (pj,Ar+i(a:,t,0), l<j<2N UjN+k = — 9^Xfi^ix,t,0), 

— 2 I 

dji = 'fij,iix,t,0), l<j<2N djk = 9 ^ 2 k^i’^ ix,t,0), 

o2fe — 2 j 

djN+i = fij^N+iixfi.O), l<j<2N djN+k = — g^l^-2*H x,t,0), 
2<k<N,l<j<2N 

The functions (p and if are defined in 

Proof : The columns of the determinants appearing in q(x^ t) are combined 
successively to eliminate in each column k (and N + k) of them the powers of e 
strictly inferior to 2(A: — 1); then each common term in numerator and denomi¬ 
nator is factorized and simplified; finally we take the limit when e goes to 0. 
Precisely, first of all, the components j of the columns 1 and -|- 1 are respec¬ 
tively equal by definition to Pji\ff\ + 0(e) for Ci, (/JjAr+i[0] -I- 0(e) for Cn+i of 
As, and + 0(e) for C[, tpjN+i[0] + 0(e) for C)v+i of Ai. 

At the first step of the reduction, we replace the columns Ct by Ct — Ci and 
CN+k by Cvr+fc — Cjv+i for 2 < A: < A, for A3; the same changes for Ai 
are done. Each component j of the column Cf of A3 can be rewritten as 

and the column Cw+fe replaced by 

l)e^* for 2 < A: < A. For Ai, we have the same reductions, each component j of 
the column C(. can be rewritten as “ 1)^^* and the column 

Cw-Hfc replaced by ■(^.i’j,N+i[l]{k‘^^ - l)e^' for 2 < A: < A. 

The term ^ for 2 < k < N can factorized in A3 and Ai in each column 
k and A -|- A: , and so these common terms can be simplified in numerator and 
denominator. 
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If we restrict the developments at order 1 in columns 2 and + 2, we get re¬ 
spectively -I- 0(e) for component j of C 2 , (pjN+i[^] + 0(e) for component 

j of Cn +2 of A 3 , and + 0(e) for component j of tpjN+ili-] + 0(e) for 

component j of C)v-i -2 This algorithm can be continued up to the columns 

Cn, C 2 N of A 3 and (7)^, of Ai. 

Then taking the limit when e tends to 0, q{x,t) can be replaced by Q{x,t) 
defined by : 


‘/5i.i[0] . 


</3i,Ar-i-i[0] 

— 1 ] 

T2,iP\ 

. ip2AN-l] 

V52,w-i-i[0] 

</22,Af-|-l[-A — 1] 

<P2Ar,l[0] . 

■ T2N,l[^ ~ 1 ] 

<P2Ar,Ar-i-l[0] . 

• T2N,N-i-l[N — M 


■ A.iW-i] 

'0l,W-l-l[O] 

— 1 ] 

"02,1 [0] 

■ tp2,l[N-l] 

'02,Af-|-l [0] 

V'2,Af-|-l[iV — 1] 

^2Af,l[0] . 

■ '4’2nA^ ~ 1 ] 

V'2Ar,Ar-i-l[0] . 

• 'lp2N,N-i-l[N — M 


So the solution of the NLS equation takes the form : 

v{x, t) = exp( 2 it — iq}) x Q{x^ t) 
So we get the result given in (Hill . □ 


4 Families of quasi-rational solutions of order N 
depending on 2N — 2 parameters 


Here a theorem which states the structure of the quasi-rational solutions to the 
NLS equation is given. It was only conjectured in preceding works [m [ 2 T]. 
Moreover we obtain here families depending on 2N — 2 parameters for the Ath- 
order Peregrine breather including families with 2 parameters constructed in 
preceding works and so we get other symmetries in these deformations than 
those were expected. 

In this section we use the notations defined in the previous sections. The func¬ 
tions ip and Ip are defined in (p^ . (l40t . (ITO . (l42l) . 


Theorem 4.1 The function v defined by 


v{x, t) = exp( 2 zt 


iip) X 


det((n^-fc),.,g,i,,„,) 


(46) 


is a quasi-rational solution of the NLS equation quotient of two polyno¬ 
mials N(x,t) and D{x,t) depending on 2N — 2 real parameters dj and bj, 
I < j < N - 1. 

N and D are polynomials of degrees N{N -I- I) in x and t. 


(45) 
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Proof : From the previous result (H5l) . we need to analyze functions (pk,i, V'fe,! 
and (pk,N+i, '>Pk,N+i- Functions (pkj and ipk,j differ only by the term of the 
argument so only the study of functions (pkj will be carried out. Then the 
study of functions ipk,j can be easily deduced from the analysis of 
The expansions of these functions in e are studied. We denote (lkj)k,je[i,2N] the 
matrix defined by 

“ Q^2j — 2 ’ ik,j+N = Q^2j — 2 ‘fk.l+N , 1 ^ J ^ N, 1 < fc < 2N, 

meaning tp. Each coefficient of the matrix ilkj)k,je[i,2N] must be evaluated, 
the power of x and t in the coefficient of column m £ [1, 27V]. We 

remark that with these notations, the matrix ilkj)k,je[i, 2 N] evaluated in e = 0 
is exactly {nkj)k,je[i, 2 N] defined in (1451) . Four cases must be studied depending 
on the parity of k. 

1. We study Iki for k odd, k = 2s + 1. 


Iki = (—1)® sin(2e(l — e^) + 4ie(l — e^) = (1 — 2e^)t 


l + 7e(l-e^) 2 , 

-iTTT-ei)xe'= "(l-e^) 

1 — 2€(1 — e^j 2 

p p 

= (- 1 )" sine(^ C2ie^‘x + 27 ^ C2ie^\l - 2e^)t + 2 ^(-l)'e' 


21 


Z =0 Z =0 

N-1 N-1 


/=0 


( 1 -e) ^ 

(27 + 1) 


- E + * E + 0(e^+')) X e'=-2(^52+"' + 0(e’-+i)) 


Z=1 Z=1 

p 


Z =1 


= (-l)®sine(^(c 2 ia; + (i 2 it + / 2 Z + 0(e^+^))e^') x e'" "(E52+"' +0(e’'+^)) 


1=0 


(_^] 1 +Sf 2 l P r 

= E (27 + 11! (E(^2-^+^2nt+/2n+0(£^+^))£"")"'+^xe^-^(^g2/£"'+0(e’^+^)) 

1^0 ^ '' n^O 1^1 

= E ^721 + 1)1 P„(x,7)e^")2'+i X e'=-i ^52+"' + 0(e‘) 

1=0 ^ h n=0 1=1 

where Pn{x, t) is a polynomial of order 1 in x and t. 


Iki = E E /3«o.....«p^o(x, tr • ■ • Ppi.X, t)«Pe2(ai+2a.+pa,) ^^2. ^ ^2+"'+0(e‘) 

/—0 CkQ-j-...-j-Ckp —2/-f-1 Z—1 


= Ee'' E xe2^^52+"' + 0(e‘), 

Z—0 cvq.- j-CKp —2Z-1-1 Z—1 

where Qao,...,aj,{x,t) is a polynomial of order 27 + 1 in x and t. 

The terms in are obtained for 7 = 0 in the two summations with ao = 1. 
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For column m, we search the terms in e2m-2 maximal power in x and 

t. It is obtained for 21 + fc — 1 = 2m — 2, which gives I = m — s — 1. 

The notations given in (|Ti)) are used. We get the following result 

Proposition 4.1 

deg(n2s+i,m) = 2(m - s) - 1 for s < m - 1, n 2 s+i,m = 0 for s > m. (47) 


2. We study hi for k even, k = 2s. 

hi = (—1)®“'"^ cos(2e(l — e^)ix + 4je(l — e^)^(l - 2e^)t 

-^7^- ^-ei) xe'^ ^(l-e^) ^ 

1 — 2€(1 — e ^) 2 

= (- 1)®+^ cose{j2 C2ie^'a; + 2t C2/e2'(l - 2e‘^)t + 2 ^ 

1^0 1^0 1^0 ^ ' 
N-1 N-1 r 

- Y, + * E + 0{eP+^)) x e'=-"(E 92 ie^' + 0(e’’+^)) 


Z=1 Z=1 

p 


Z=1 


= (-l)'*+^cose(^(c2za; + d2it + /2i + 0(e^’+^))e^') x e’" g2 ie^'+ 0(e’'+^)) 


1=0 


1=1 


— E -Wv-(E(^2na^+c^2nt+/2n + 0(e^''’^))e^"')^^xe^ ^(E 92lf^''+Oh^~^^)) 


( 21 )! 

l—O ^ ' n—0 


Z=1 


= E ‘ (E «)*"")” « + o('‘) 


( 21)1 

Z=0 ^ n=0 

where Pn{x^ t) is a polynomial of order 1 in x and t. 


z=i 


^^4 = E E Po.o,-,c.,Po{x, tr • . • Pp{x, t)“Pe2(a,+2a.+pa,) ^^2.-2 ^ 

Z—0 ckQ-j-.. .-j-ctp—2Z Z=1 


= E E (^, t)e2(“i+2“^+P“r) x £2-2 ^ ^ 

Z—0 ctQ “1“ ■ ■ ■ “l“Q^p —2Z Z —1 

where Qao,...,ap{x,t) is a polynomial of order 21 in x and t. 

The terms in e° are obtained for 1 = 0 in the two summations with oq = 1. 

For column m, we search the terms in e 2 m -2 .y^rith the maximal power in x and 
t. It is obtained for 21 + fc — 2 = 2m — 2, which gives I = m — s. 

With the notations given in (l44l) . we have 

Proposition 4.2 

deg(n 2 s,m) = 2(m - s) for s < m, n^s^ra = 0 for s > m. (48) 
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3. We study lkM.+i for k odd, k = 2s + 1. 


= (—1)® cos(2e(l—e^) 2x—4ie(l—e^) 2 (1—2e^)t+i In 


l + ze(l-e^) 2 


-5-CM , -I / 

l-ze(l-e 2 )-^ 


xe"-'=-i(l - '2^-- 


= (—l)®(cose(^ C2/e^^a; — 2z ^ C2;e^^(l — 2e^)t — 2^(—^ ^ 


2^-- 


1=0 1=0 
N-l N-1 


1=0 


{21 + 1 ) 


- ^ + z ^ bie^^ + 0(eP+i)) x g^-fc-i ^ o{e^+^)) 

1^1 1^1 1^1 

p r 

= {-lyicos t{Y,{c 2 ix + d 2 it + f 2 i)e^‘ + 0{eP+y) x g 2 ie^‘ + 0{e^+y) 


1=0 
]!+s^2l P 


1=1 


E ^ (E(^2„x+d2„t+/2„+0(e^+^))e^")^^xe^-^-^(^g2K^^+Q(e-+^)) 


/=0 

q 


n—O 

^21 P 


/=1 


= E MS^(E + 0{eP+yr X 


.M-2S-2 


(20! 

/=0 ^ ^ n=0 


iJ292ie^' + Oie^+y) 


1^1 


where Pn(^5 0 ^ polynomial of order 1 in x and t. 


q 

lkM + i=J2^^^ E Pcco,...,a^Po{x,tY° 

/—0 ocQ-\-.. .-\-Oip—2l 


Pp{x i)“Pg2(ai+2a2+pap) 


xe'^-2--2^52ie'' + 0(e‘) 


Z—0 ctQ + ,, ,+ctp— 2/ 


xe'^-2^-2^52ie2' + 0(e‘), 
/=! 


where Qao,...,ap{x,t) is a polynomial of order 21 in x and t. 

The terms in (column -^ + 1) are obtained for ^ = 0 in the two summations 
with ao = 1- 

For column ^ + m, we search the terms in e2iTi-2 maximal power in x 

and t. It is obtained for 21 + 2{N — s—1) = 2m — 2, which gives I = m + s — N. 
Then we get the following result 


Proposition 4.3 

/ N M M ^ ^ 

d(ig{n 2 s+i,m+Jf) =‘2m + 2s- M for s> — - m, n 2 s+i,m = 0 for s < — - m. (49) 
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4. We study lk,i+^ ^ even, k = 2s. 

= (— 1 )^ sin( 2 e(l—e^)^x—4te(l—£^)^( 1 — 2 e^)t+i In ^ 


1 — je(l — e 


2'>-4 T 


xe 


M-k-l 


, On M -K- 

(l_e2)-r- 


P P P n - ^2',-^ 

= (-1)® sine(^ C2/e^'a; - 2i'^C2ie^\l - 2e^)t - 2^(-l)'e^' 
i=0 1=0 1=0 


{21 + 1 ) 


N-l 


N-1 


- Y, + 0{eP+^)) X + 0(e'-+i)) 

p r 

= (- 1 )® sin e(^(c 2 /:r + d 2 it + f 2 i)e^^ + 0{eP+^)) x ^M-fc-i ^ o{e'^+^)) 

1^0 1^1 

^ P 'T 

- E Vn/ , 1 V (E(^2„a^+d2nt+/2„+0(£"+^))£^")^^+^xe^-n^g2/e^^+0(e-+^)) 


(2Z + 1)! 

1=0 ' ^ n=0 

(- 1 )' 


1=1 


= E E+li! (E X e^^-^®(^g2/e^' + 0{e^+^)) 

1=0 ' '' n=0 1=1 

where Pn{x,t) is a polynomial of order 1 in a; and t. 

9 

lk.i=Y^^’’ E /3ao.....ap^o(a:,t)““ 

l—O CtO“l“---“|“Q;p — 2 /- 1-1 

r 

. . . Pp{x, t)“Pe 2 («i+ 2 a 2 +pap) ^ ^M- 2 s ^ 

1^1 

= E^'' E Qco....,a,(x,t)e 2 (“i+ 2 “^+P“r) 

/—0 olq-\-.. .-\-(y.p^‘2,l~\-\. / — I 

where Qao,...,ap{x,t) is a polynomial of order 2 ^ + 1 in a: and t. 

The terms in are obtained for / = 0 in the two summations with = 1. 

For column ^ + m, we search the terms in e 2 m -2 maximal power in x 

and t. It is obtained for 21 + M — k = 2m — 2, which gives l = m + s — A^ — 1. 
Using the notations given in (H^ . we get the following result 


Proposition 4.4 


deg(n 2 ^ m+M.) = 2m + 2s — M — 1 for s > — + 1 — M, 

^ M 

2 


^ 2 s,m+ 


M = 0 for s < ^ -\- 1 — m. 


(50) 


These results can be rewritten in the following way 
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Proposition 4.5 


deg(nj,fe) = 2k-j for j < 2k, 

Uj^k = 0 for j > 2k, , , 

deg{nj^k) = 2k + j — 2M — 1 for j > 2M + 1 — 2k, ' 

= 0 for j < 2M + 1 — 2k. 

The degree of the determinant of the matrix {nkj)k,je[i, 2 N] can now be evalu¬ 
ated. 

From the previous analysis, we see that x and t have necessarily the same power 
in each rikj- The maximal power in x and t, is successively taken in each column. 
It is realized by the following product 

N N 

n n nN+j, 2 N+i-i- 

t=i i=i 

Applying the result given in ([ST|) we get 

N N 

deg{det{nkj)k,j(^li,2N]) = deg(njj) -b deg{nN+j,2N+i-j) 

t=i j=i 

N N 

= ^ 2j - J + ^ 2(M + 1 - j) - 2M - 1 + — + J 
i=i i=i 

N N 

= ^j + ^iV+l-j = iV(iV+l). 
t=i i=i 

It is the same for determinant det{dkj)k,j&[i,2N], we have deg(det((i/cj)j, 2Af]) = 
N{N + 1). 

Thus the quotient 

det((<ifej)^. jAT]) 

defines a quotient of two polynomials, each of them of degree N{N + 1), and 
this proves the result. 

Parameters oi = J2k=i dfcCfe ^md oi = X]fc=i^ dkCk must be chosen in the fol¬ 
lowing way. 

The term Ck must be a power of e to get a nontrivial solution; Ck must be a 
strictly positive number a in order to have a finite limit when e goes to 0. If 
the power of e is superior to 2N — 2, the derivations going up to 2N — 2, then 
this coefficient becomes 0 when the limit is taken when e goes to 0 and so has 
no relevance in the expression of the limit. 

□ 
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5 Conclusion 


Here we proved the structure of quasi-rational solutions to the one dimensional 
focusing NLS equation at order N. They can be expressed as a product of an 
exponential depending on t by a ratio of two polynomials of degree N{N + 1) in 
X and t. If we choose = 0 for 1 < z < iV — 1, we obtain the classical (ana¬ 

logue) Peregrine breather. Thus these solutions appear as 2N — 2-parameters 
deformations of the Peregrine breather of order N. 

The solutions for orders 3 and 4 first found by Matveev have also been explic¬ 
itly found by the present author [26j W\- We have also explicitly found the 
solutions at order 5 with 8 parameters [28]: these expressions are too extensive 
to be presented : it takes 14049 pages! For other orders 6, 7, 8, the solutions 
are also explicitly found but are too long to be published in any review. In the 
relative works dB [13 uni inni Em only the analysis has been done and figures 
of deformations of the Peregrine breathers has been realized. The solutions for 
order 9 with 16 parameters [33] and respectively for order 10 with 18 parameters 
are also completely found [34] . 

We still insist on the fact that quasi rational solutions of NLS equation can 
be expressed as a quotient of two polynomials of degree N{N -|- 1) in a: and t 
dependent on 2N — 2 real parameters by an exponential depending on time. 
Among these aforementioned solutions of order N, there is one which has the 
largest module : it is the solution obtained in this representation when all the 
parameters are equal to 0; one obtains the Peregrine breather order N. His 
importance is due to the fact that among the solutions of order iV, its module 
is largest, equal to 2N + 1. This result first formulated by Akhmediev has just 
been proved recently [55] . 

In the recent studies proposed by the author, the solutions of order N can be 
represented by their module in the plane (a;;t). With the representation given 
in this article, one obtains at order N, the configurations containing N{N -|-1)/2 
peaks, except the special case of Peregrine breather. These configurations can 
be classified according to the values of the parameters Oi or bi for i varying be¬ 
tween 1 and — 1. It is important to note that the role played by Oi or bi for a 
given index i is the same one, in obtaining the configurations. The study refers 
to the analysis of the solutions when only one of the parameters is not equal to 
0. Among these solutions, one distinguishes two types of configurations; for oi 
or bi not equal to 0, one observes triangular configurations with N{N + l)/2 
peaks. For or bi not equal to 0 and 2 < * < — 1, one observes concentric 

rings. The simplest structure is obtained for ojv-i or 6jv-i not equal to 0 : one 
obtains only one ring of 2N — 1 peaks with in his center Peregrine breather of 
order N — 2; this fact was also first formulated by Akhmediev. The detailed 
study of the other structures is being analyzed. We hope to be able to give 
results soon. 

We can conclude that the method described in the present paper provides a very 
efficient and powerful tool to get explicit solutions to the NLS equation and to 
understand the behavior of rogue waves. 

There are currently many applications in different fields as recent works by 
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Akhmediev et al. [35] or Kibler et al. m attest it in particular. 

This study leads to a better understanding of the phenomenon of rogue waves 
and it would be relevant to go on with higher orders. 
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